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abstract 

A joinc experiment becween Canadian and USA 
research teams was conducted early in October, 
1992 to determine the effect of species structure 
and dielectric variations on forest backscatter. 

Two stands, one red pine and one jack pine. «■« c e 
Petawawa National Forestry Institute (PNFI) were 
utilized for the experiment. Extensive tree 
architecture measurements had been taken by e 
Canada Centre for Remote Sensing (OCRS) several 
months earlier by employing a Total Station 
surveying instrument which provides detailed 
information on branch structure. A second part of 
the experiment consisted of cutting down several 
trees and using dielectric probes to measure 
branch and needle permittivity values at both 

S1CeS The dielectric and the tree geometry data 
have been used in the Ceorge Washington University 
(CUU) Vegetation Model to determine the C band 
backscattering coefficients of the individual 
stands for W polarization. The model results show 
chat backscatter at C band comes mainly from the 
needles and small branches and che upper portion 
of che trunks acts only as an attenuator. The 
paper will conclude with a discussion of variation 
of backscaccer with specie structure and how 
dielectric variations in needles for both species 
mav affect che total backscatter returns. 
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INTRODUCTION 

A joint experiment between Canadian and USA 
research teams was conducted late in September, 
1992 to determine the effect of species structure 
and dielectric constant on forest 
backscatter. Two stands, one red pine and one jack 
pine, in the Petawawa National Forestry Institute 
(PNFI) were utilized for the experiment. C and X 
band dacs was obtained from over-flights by the 
OCRS airborne SAR and measurements were made of 
tree architecture and branch and needle dielectric 
constancs. This daca was then used in the CUU 
vegetation model so that model and measurement 

results could be compered. 

The two test sites were located Inside the 
PNFI research forest which is centered at <*6 00‘N 
77*25 ‘W, near Chalk River, Ontario, Canada. Both 
test sites were characteristic of plantation 
stands having even age trees which were planted on 
flat ground in a regular fashion. The red pine 
stand CPinus resinosa) , which has an extent of . 
hectares, consisted of 58 yr old tr..s with an 
average diameter at breast height (dbh) of 23_5 
and a 8 height of 21 m. The jack pin. stand (Pinus 
ban Asians) consisted of 26 years old trees which 
had an average dbh of 20 cm and a height of 15 m, 


they covered a 3.5 hectare area. The previous July 
(1992). as part of an optical experiment, care.u 
architectural measurements using a vectorizatton 
technique were made of three trees in each stand. 

At the time of the over-flights dielectric 
measurements of branches and needles were made at • 

C and X bands by using Applied Microwave porcable 
dielectric meter with tree probes. The CuRi 
airborne SAR flew over the sites on September 25 
and 26 1992 The flights were coordinated witn 

calibration test flights OCRS was performing and 
thus comer reflectors and active radar 
calibrators (ARCs) were in place for calibration 
purposes. On September 25th Like and cross 
polarized data was obtained at C (5.3 GHz) and X 
(9.25 GHz) band for several angles of Incidence. 

On the following day, polarimetric data was 
collected at C band. The weather at the time or 
both over-flights was clear and dry In the 
present paper only C band W polarized data will 

be employed. , . 

Several researchers have been interested 
tree structure and dielectric constant as a 
mechanism for using radar to classify trees y 
specie. In the literature there l* some evidence 
that such differences exist. Leckie [l], using CCRS 
C and X band SAR data has demonstrated that 
measurable differences exist between backscatter 
values for different tree types Inthe ™FI 
forest. Dobson, et al(2] . employing jPL SAR P, L 
and C band measuremants made in North Michigan, 
has shown that both deciduous and conifer stands 
have discernible backscatter differences over the 
spring to summer time-frame. Lang et al[3] has 
shown how the long straight trunks of the red pine 
are partially responsible for the bright ” CU ™ 
observed at P band. Salas, et al(4] and McDonald 
et all 5] have demonstrated how the dielectric 
constanc in tha accive xylem responses to water 
conditions in the tree. These changes in 

dielectric constanc in che outer sapwood may 
effect backscatter values particularly ac C and X 
bands whara the penetration of electromagnetic 
waves is small. 

tree architecture measurements 

Tree architecture data was collected using a 
new technique called ’Tree Vectorization 
developed ac the Canada Cencre for Remote Sensing^ 
The method was developed to characterize i 
architecture from a restricted daca tofU ' 
use this sample to reconstruct a statist! 
accurata 3 -dlmensionaL representation ofch. ere 
Landry et al [6). The methodology can be use 
for cree and forest canopy modeling research l 
both microwave and optical part o 
spectrum. 



The vector izat ion method involves data 
collection at three different levels: the trunk 
inventory, the branching structure, and the 
foliage data collection, each has a specific 
methodology and provides a separate set of data. 
The spatial distribution of the supporting 
r. r rhe primary branches, is characterized through 
the trunk inventory The branch structure is 
specified by a vectorized sampling method, and 
concurrently, foliage data is recorded for 
individual branch components. Finally, the 
information, which is collected at all levels, is 
combined at the final stage and used for tree 
reconstruction. 

Briefly, the sampling strategy proceeds by 
first making a trunk inventory. The spatial 
distribution of all major branches that leave the 
trunk along with their diameter and elevation 
angle are recorded. The branch structure is 
defined by sampling a sufficient number of 
primary branches to get a representative variation 
of the branch structure along the tree height 
axis . 

The structure of these selected branches is 
characterized in terms of branch segments. Each 
segment corresponds to that portion of the branch 
between two branch nodes or ramification points. A 
bclecLed number of connected branch segments or 
paths are designated for sampling and the spatial 
location of each segment's nodes is recorded using 
a survey instrument (Vild TC1000 Total Station). 
The precision of this instrument is more than 
adequate for these purposes since it is within 3 
ram in a 2 km range, and within 3 seconds of arc in 
horizontal and vertical angles. Each sampled 
segment is tagged and attributes such as segment 
number, connecting segment numbers, diameter and 
foliage information are recorded. 

From the sampled branch information and the 
trunk inventory, the 3 dimensional architecture of 
the complete tree can be reconstructed. The 
reconstruction is not exact since similarity 
principles are employed to reconstruct whole 
branches from the sampled portions. The 
reconstructed branches are then used to represent 
other primary branches on the trunk that were not 
selected for sampling. In this way an approximate 
representation of the tree in terms of segments 
whose coordinates are known is obtained. The 
method is discussed in more detail in Landry [6], 

The vectorized sampling and reconstruction 
technique has been used to model one typical tree 
in the red pine stand and one in the jack pine 
stand. From the reconstructed trees, segment data 
has been used to determine branch paths. Since the 
segments In each branch are not necessarily 
colllnear, the reconstructed branches are, in 
general, curved. Since branches used in the 
scattering model are cylinders, the reconstructed 
branches were replaced by cylinders that have the 
mean diameter and the mean elevation angle of the 
reconstructed branch. These averages are defined 
by weighting each segment’s diameter or elevation 
angle by its length These mean branches or 
cylinders are then classified into size categories 
based on their diameter and branch length 
distributions ; these distributions are determined 
from all the recons tructed branches in the tree. 
The branch classes are listed in Table 2. A twig 
class has been created since twigs don’t usually 
occur at nodes and thus must be treated specially. 
Also, the branch class 0 tends to reflect the 
terminal segments found in the tree. The terminal 
segment class and twig class have similar 
statiscicc in terms of mean diameter and length. 

The probability density functions of branch 


length. mean branch diameter and mean branch 
elevation angle and the density of branches per 
cubic meter were derived for each branch class 
The number of needles was calculated by using the 
information collected per segment TtocaL length of 
segment with needles, percentage of needles and 
number of r.czdles pec unit u! length;, ana tnen, 
by summing segment values over the whole tree. 

DIELECTRIC CONSTANT MEASUREMENTS 

Detailed measurements of the dielectric 
constant for a red pine and a jack pine tree were 
made using a dielectric probe. The use of a 
dielectric probe to obtain dielectric constant 
values was preferred over the empirical approach 
used by the authors in the past since the probe 
measures the dielectric constant directly. The 
empirical approach estimates the dielectric 
constant and this estimate is primarily dependent 
upon moisture concents of the vegetation. 
Important electromagnetic parameters, such as 
ionic conductivity is usually assumed to be the 
same for all parts of a tree. The aim of the 
probe measurements, in the present experiment, was 
to investigate: 1) differences in the dielectric 

constants of different species, 2) differences in 
the dielectric constant in the same species. At 
the time of the experiment, it was thought that at 
C band, the trunk of a tree is less important chan 
the crown. Therefore, the dielectric constant 
measurements were limited to the branches and the 
needles only. The measurements on branches were 
done in the fields, whereas, the clusters of 
needles were sealed in plastic bags and their 
dielectric constant was measured in the 
laboratory. 

The live crown of a red pine and a jack pine 
tree was divided into three segments. A typical 
branch from each segment was selected. The branch 
was cut into three parts. To study the 

differences in dielectric constant as a function 
of depth inside a branch, the cross-section of a 
branch was divided into three annular regions. 
The region just beneath the bark was called outer 
region. The center part, looking into the 
cross-section of branch was called inner region. 
The remaining area between the outer and inner 
regions was called middle region. Probe 

measurements across the cross-sec tion of a branch 
segment were taken in each of the regions. It was 
found that values for the center and middle 
regions were of the same order. As a result, the 
center region was expanded to encompass the middle 
region to simplify the discussion. Fig. 1 shows a 
typical branch cross section with the two regions 
labeled. The outer region has a width L^ while 

the radius of the center region is L^ . 



Dielectric measurements over the 
cross-section cf a branch 



As an illustration, the real part of the 
measured dielectric values from a red pine tree 
for the cencer and outer regions are plotted in 
Fig 2 The imaginary part is found to follow the 
same trend, and therefore is not shown in the 

figure. The differences in the dielectric constant 
(both real and imaginary pa- -s ) ~ 

and cencer regions exist for branches of all 
diameters. The high values of dielectric constant 
that are observed in the outer region are most 
likely due to sap and water transport layers chat 
exist there. The region has a thickness, , of 


3mm This thickness is determined by the probe 
tip width. It should be noted from Fig. 2 that the 
outer and center dielectric appear to converge as 
the diameter of the branch becomes smaller. This 
is because che minimum resolution of the probe tip 
makes it increasingly difficult to resolve outer 
and cencer regions for small branch diameters. It 
should be mentioned here that the data from all 
segments of the tree are merged together to obtain 



0.5 1 1.5 2 2 5 3 3 5 4 

0R*CM OWA46TER 04 

Fig. 2 

Dielectric variations with branch diameter 


Ir\ order to measure the dielectric constant 
of needles, a hose clamp was used to bundle the 
needles tightly together. The clamp was squeezed 
so as to remove the air gaps among the needles . 
Care was taken not to crush the needles by 
over-squeezing the damp. The bundle was cut in 
half to obtain a flat face of the needles. A mean 
of five probe readings taken across the 
cross-section of the bundle was used to specify 
the dielectric constant of the needles. Four 
bundles of needles were used for dielectric 
constant measurements. 


average results from branches ar.d needles for red 
pine and jack pine trees. The species difference 
are clearly noticed in the case of branches, 
while, for needles such differences are not found. 
We also determined wet and dry weights of the 
branch segments and the needles The moisture 
contents were determined from this data for both 
species. No significant differences in moisture 
contents were discovered between che red pine and 
jack pine '"aes. 

FOREST MODEL 


The forest is modeled as a two layered random 
medium with a rough surface beneath. The tree 
components such as branches, needles, and crunks 
are distributed in the medium according to their 
measured orientation statistics. The distorted 
Born Approximation is used to compute the total 
backscatter lng coefficient, a , for the forest. 

The total return can be expressed as a sum of the 
direct, direct-reflected and surface backscatter 
terms; however at C band the main contribution is 
observed to be from the direct term only. The 
geometric parameters of the forest medium; such as 
density, scatterer sizes ’and orientations are 
obtained from the tree veccorizat ion data. For 
modeling puLposes, che small size branches and che 
twigs of the vectorization method were classified 
as secondary, and che branches which were larger 
were named as primary. The probability density 
function for che needles for both red pine and 
jack pine stands were chosen as sin 23, for 0<tf<90 
degrees . The secondary and primary branches for 
these stands were assigned uniform distributions 
whose end points were estimated from the 
vectorized data. 

The dielectric measurements were made for the 
outer and center parts of the branches , as seen in 
Fig. 1. The data for these regions was used to 
determine che dielectric constant of an 
equivalent, homogeneous branch of the same size. 
As seen in Fig. 1. the outer part of che branch, 
which has the highest dielectric value is assumed 
to be as thick as the probe tip. The thickness of 
the center and outer parts are denoted by . and 

L respectively. t and i represent che 

out J ct out 

dielectric constants for these regions and are 
determined from the measurements . An average 
dielectric value, t for an equivalent branch of 

• V 

che same size is then calculated by using 


k n L - k n L 

o AT O out out 


+ k n L 

o ct ct 


( 1 ) 


TABLE 1 

Comparison of Mean Values of Real and Imaginary 
Part of Dielectric Constant 


Scatterer Type 

Red Pine 

Jack Pine 


f Center 

(16.7, 

3.9) 

(3.8, 2.0) 

Branch 

1 

[ Outermost 

(37.1, 

16.7) 

(25.3, 10.5) 

Needles 


(13.4, 

5.4) 

(13.3, 5.6) 


It took 4 days to complete dielectric 
constant measurements. The measurements for 
branches of a red pine tree were made mostly 
during the mornings and the jack pine data were 
collected In the afternoon. The calibration 
accuracy of the probe was checked periodically by 
using standard liquids such as butanol, ethanol, 
methanol, formaaido, etc. Table 1 summarizes 


where n is the refractive index and is given by 
the square root of the dielectric constant. The 
right hand side of Eq. (1) represents the complex 
phase from the branch surface to its center. This 
is equated to the complex phase for a homogeneous 
branch of permittivity It should be noted 

that this approach will be valid only for the 
cases when the skin depth exceeds L. As the skin 
depth gets smaller, the dielectric constant should 
approach the dielectric value of che outer part. 
However, Eq . (1) will not differentiate between 
these, cases and will underestimate the average 
permittivity when the skin depth is less than L. 
Table 2 below summarizes the values for the model 
parameters for red pine and jack pine stands as 
derived from the vectorized geometric data and Eq . 
(1). The trunk layer has been excluded since at 
C band backscattering coefficient originates 
mainly from che tree crown. 




TASLS 2 

Modal Partaacan for Rad Pina and Jack Ptna Traaa 
Rad Pina 


Scec. Kind 

Length 

Rad. 

Oanalcy 


Ty?« 

'em) 

(ca) 



NEEDLE 

10.30 

0.06 

7786.0 

3 4+L5 . 4 

5tC.BR. cvl* 

4 20 

0 32 

3 13 

6.7+115. 3 

br. 0 

19.30 

0.33 

6.13 

6.7+115. 3 

br . 1 

26.50 

0.47 

9.48 

6,7+112.2 

?RI BR. br. 2 

105.50 

0.75 

0.30 

6.2+17.3 

br . 3 

2 34.00 

1.05 

0.38 

13. 5+16. 6 


Jack Pine 



Scac. Kind 

Langch 

RAd. 

Density 

< 

A* 

Type 

(cm) 

(ca) 

<•/■’) 


NEEDLE 

3.10 

0.05 

39360.0 

13.4+15.6 

SEC.BR. cvlg 

7.35 

0.14 

50.9 

25.0+110.0 

br . 0 

8.29 

0.L0 

91.80 

25.0+110.0 

br . 1 

19.00 

0.14 

70.10 

25.0+110.0 

PRI. ML br. 2 

32.10 

0.24 

27.50 

25.0+110.0 

br . 3 

53.00 

0.31 

2.11 

20 . 0+18 . 0 

br. 4 

172.50 

0.68 

1.55 

15.0+16.0 


RESULTS AND DISCUSSION 


W 


che 


In chis section che attenuation and 
backscatter predictions of che GVU Vegetation 
Model for both red pine and jack pine stands will 
be examined. Also. che reasons for different 
backscatter returns for these species and their 
sensitivity to dielectric variations will be 
analyzed. The discussions will involve 
polarization only. 

Using che parameters given in Table L, 
attenuation due to each scatterer type in both red 
pine and jack pine stands have been estimated. As 
it can be seen from Table 3. che main factor 
causing attenuation in the red pine forest is the 
needles; whereas for the jack pine both needles 
and the primary branches play an equally important 
role . 

TABLE 3 

Total Incoming Attenuation (dB) 
for W Polarization 

Scat. Type EiM. 


NEEDLE 
SEC. BR. 
PRI. BR. 
Total 


5.68 

3.48 

5.73 

14.89 


6.90 

3.00 

1.46 

11.36 


The model predictions for the backscatter 
coefficients for each of these species is given in 
Table 4. If the individual returns for each of the 
scatterer components are compared, one can observe 
that the major contribution in the red pine stand 
is obtained from the needles; whereas for the jack 
pine forest it is the secondary branches that give 
the highest returns. The fact that the needles 
have che highest return in the red pine stand can 
be explained by the low density of the secondary 
branches, and also by long needles. On che other 
hand; che density of the secondary branches in 
jack pine is much higher compared to the red pine 
and thus their contribution to backscattering 
coefficient can surpass chat of needles which are 
very small in size. 


The sensitivity of the backscatter returns to 
che dielectric constant of needles is examined in 
Fig. 3. One can observe the variations in the 
total backscattering coefficient as the dielectric 
constant of the needles for both red pine and jack 
pine stands is varied from %50 to %300 of their 
nominal values, given in Table 2. It is seen that 
the red pine returns increase by as much as 5 dB , 
while che response from the jack pine stays almost 


constant . 


TABLE 4 

Backscatter Coefficient (dB) for W Polarization 



Kind 

Red Pine 

Jack Pine 

NEEDLE 


-17.99 

-22.53 

SEC. BR. 

tvig 

-32.93 , 

_ -i ■) ^ 


hr . 0 

-25.98 

-22.86 


br . 1 

-26.63 ! 

-15.32 

PRI. BR. 

br. 2 

-36.76 

-18 . 77 


br . 3 

-20.95 

-38 . 54 


br. 4 

j 

-20.04 

Total 


-15.33 

-LI. 68 

Measured 


-13.55 

| -12.07 



Ratio of the change in the dielectric constant 
from che nominal value (%) 
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